Proposed search for T-odd, P-even interactions in spectra of chaotic atoms 
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Violation of fundamental symmetries in atoms is the subject of intense experimental and theo- 
retical interest. P-odd, T-even transitions have been observed and are in excellent agreement with 
electroweak theory. Searches for permanent electric dipole moments have placed bounds on T-odd, 
P-odd interactions, constraining proposed extensions to the Standard Model of elementary parti- 
cles. Here we propose a new search for T-odd, P-even (TOPE) interactions in atoms. We consider 
open-shell atoms, such as the rare-earth atoms, which have dense, chaotic excitation spectra with 
strong level repulsion. The strength of the level repulsion depends on the underlying symmetries 
of the atomic Hamiltonian. TOPE interactions lead to enhanced level repulsion. We demonstrate 
how a statistical analysis of many chaotic spectra can determine the strength of level repulsion; in 
particular, the variance of the number of levels in an energy range has been shown to be a useful 
measure. We estimate that, using frequency comb spectroscopy, a sufficient number of chaotic levels 
could be measured to match or exceed the current experimental bounds on TOPE interactions. 
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I. BACKGROUND 

Electromagnetic interactions are invariant under any 
combination of the discrete symmetry operations, includ- 
ing spatial parity (P), time- reversal (T), and charge con- 
jugation (C). Conservation of these symmetries lead to 
the well-known selection rules for atomic expectation val- 
ues and transition amplitudes. Searching for violation 
of these selection rules enables revealing minute non- 
electromagnetic corrections to the atomic Hamiltonian. 
First observation of the parity-forbidden El transition in 
bismuth was an important confirmation of the electro- 
weak theory in semi-leptonic sector Later atomic 
experiments with cesium confirmed low-energy limit of 
the electro-weak theory at the level of radiative cor- 
rections 043|- While the parity violation is firmly es- 
tablished, the T-violation is probed in searches for the 
electric dipole moments (EDMs) of elementary particles, 
which violate both P and T symmetries @. The most 
stringent limits to date on the electron EDM were derived 
from experiments with atomic thallium Q and with YbF 
molecule Q. The most accurate limit on atomic EDM 
comes from Hg experiment [8]. It is worth emphasizing, 
that while no EDM has been discovered yet, the current 
limits on eEDM are very powerful. These limits have 
ruled out many extensions to the Standard model in- 
cluding minimal supersymmetry and severely constrained 
SUSY models in general. 

Here we focus on the T-odd, P-even (TOPE) interac- 
tions. The TOPE interactions are far less explored than 
the described P-odd and P,T-odd couplings. One rea- 
son for this paucity is of theoretical nature: phenomeno- 
logical Lagrangians for TOPE interactions involve field 
derivatives; thereby in most popular models such terms 
appear only as radiative corrections. Another reason 
is purely experimental: there is no convenient observ- 
able, like EDM, associated with TOPE interactions. This 
makes designing experiments more challenging. So far 
only two types of atomic physics experiments were sug- 



gested and carried out for TOPE interaction searches. 
Here we propose an alternative to these searches. 

Before describing our proposal, we review the previous 
efforts and the established limits. The positronium ex- 
periment Q relied on the compelling arguments of the 
CPT theorem which implies that non-vanishing TOPE 
interactions would violate charge conjugation symmetry. 
The experiment probed the C-forbidden single-photon 
transition between the 2 3 Si and 2 1 Pi states of positro- 
nium. That lead to a constraint on the electron-positron 
TOPE interaction H°~ B : 
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The second search [l0( , proposed by Kozlov and Porsev 
ll| , focused on the T-odd correlation k ■ E in the refrac- 
tive index of atomic vapour near the 6pi/2 — > &P3/2 
transition in thallium. Here k is photon momentum and 
E is external DC electric field. The experiment on Tl 
was carried out on a particular hyperfine transition and 
placed a limit on the electron-proton interaction (with 
valence proton in Tl nucleus). One can think of this in- 
teraction as a T-odd hyperfine structure. 

The experiment [Ic| placed the following limit on the 
nuclear-spin-dependent (NSD) interaction: 



x 1! - T ' '-' = (0.9 ± 2.0) x 10" 



(2) 



Up to now ([T]) and @ are the only direct experimental 
limits on TOPE interaction in atoms. Stronger indirect 
model-dependent limits can be obtained from the EDM 
experiments. Indeed, in the second order of perturbation 
theory the T-odd and P-odd interactions generate effec- 
tive P,T-odd interaction. This, in turn, contributes to the 
EDMs of elementary particles and atoms. Khriplovich 
noted that for the short-range interaction this mecha- 
nism leads to very stringent bounds on the short-range 
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TOPE interaction [I1GJ]. For the long-range TOPE in- 
teraction such indirect limits are less stringent, but still 
stronger, than the direct limits (11121) (Toj |. Note, how- 
ever, that the indirect limits are mo del- dependent, i.e., 
they depend on assumptions with respect to employed 
phenomenological Lagrangians. 

In nuclear physics, the constraint on TOPE comes 
from level-repulsion statistics in chaotic nuclear spec- 
tra [3, [HI . For such a system, the Hamiltonian may be 
represented as a block-diagonal matrix, where each block 
corresponds to states of the same total angular momen- 
tum and parity. If we are merely interested in statistics 
of level spacings, we can avoid diagonalizing or even con- 
structing the Hamiltonian. Instead one could consider 
ensembles of random matrices whose matrix elements are 
independent random variables (for nuclear physics appli- 
cations, see pjj]). These random variables are described 
by Gaussian distribution with a fixed variance. Ensemble 
averages of fluctuation measures (such as the E 2 measure 
used in Ref. [lij]) agree well with experimental data, even 
despite unrealistic many-body interactions such matrices 
contain. 

More generally, any open-shell system of interacting 
fermions exhibits quantum chaos. Technically, any quan- 
tum system with a nonseparable wave equation is chaotic; 
this includes even such simple systems as neutral he- 
lium [l7| ■ For our purposes, this is much too broad a 
criteria. We are interested in atoms with chaotic spec- 
tra that are also dense (containing a statistically large 
number of levels). We also require strong interactions 
between valence electrons so that total angular momen- 
tum and parity are the only "good" quantum numbers 
(see Sec. IIIII as well as [l8| for further discussion of this 
point). 

For these reasons, the rare-earth atoms are the ideal 
candidates for our search. Their Hamiltonians can be 
modeled with random matrix ensembles exactly as de- 
scribed above for nuclei. Practically, the most useful 
criteria to distinguish between chaotic and regular (in- 
tegrable) quantum systems is the distribution of near- 
est neighbor spacings of energy levels. Regular systems 
exhibit an exponential distribution for nearest-neighbor 
spacings, while chaotic systems follow the Wigner distri- 
bution [171. 1 1 81] . 

The basic idea of this paper is to point out the utility 
of the level-repulsion statistics method of nuclear physics 
for searching for TOPE interactions in atomic physics. 



II. METHOD 

Spectra in nuclei and complex open-shell atoms exhibit 
quantum chaos, characterized as described in Sec. |U We 
can model such many-particle systems with random ma- 
trix ensembles (RMEs). If T- reversal is preserved, all the 
Hamiltonian matrix elements may be taken as real ran- 
dom variables, whereas if T-reversal is broken, the matrix 
elements are in general complex. The real and imaginary 



parts of Uij are independent random variables. There- 
fore, as described in the appendix, T-violating interac- 
tions lead to enhanced level repulsion (quadratic, rather 
than linear, for small spacings). 

To search for TOPE interactions, following [3], we use 
a RME of the form 

{H a } = {H(S)+iaH{A)}. (3) 

Here the {H(S)} matrix is real and symmetric, repre- 
senting the dominant electromagnetic interactions, and 
{H(A)} is real and anti-symmetric, characterizing TOPE 
interactions. We assume that both random matrices have 
the same variance v 2 for all matrix elements, with the 
real parameter aCl determining the relative strength 
of TOPE and electromagnetic interactions. The a = 
case is conventionally termed the Gaussian orthogonal 
ensemble (GOE), while a = 1 is known as the Gaussian 
unitary ensemble (GUE), in reference to the symmetry 
groups their respective Hamiltonians posess. From per- 
turbation theory, it should appear plausible that A, de- 
fined by 

A = a 2 v 2 /D 2 {E) (4) 

is a more useful parametrization of the relative strength 
of the TOPE interactions than a; one may think of it 
as the square of the expected value of the perturbation 
(although technically, (H(A)) = since H(A) is anti- 
symmetric). 

Before we can compare experimental spectra with 
RME predictions, we must "unfold" the spectra. This is 
necessary because RMEs have a uniform level density (al- 
ternately, the average level spacing is uniform) whereas 
real spectra clearly do not. To unfold the spectra, one 
would first construct the number staircase function A/". 
Formally, it is given by [16| 

N{E)= j E ^2 I 5{E' — Ei) dE', (5) 

i.e., M(E) is simply the number of levels with energy less 
than or equal to E. By fitting a polynomial of reasonable 
order to M(E) : one would obtain Nfu{E) : and differenti- 
ating this with respect to E would give the average level 
density p(E) [32l |. The final step would be to map each 
old eigenvalue Ei to a new eigenvalue with the for- 
mula Si = J E ^p(E)dE. But since p(E) = -^N flt {E), 
the unfolded eigenenergies Ei are simply given by 

e i =Aff it (E i ), (6) 

so that the unfolded energies are dimensionless with uni- 
form average level spacing equal to unity. 

We now introduce statistical measures with which we 
may compare unfolded data to the predictions of the 
RMEs. For our purposes, the most useful statistic is 
S 2 (r), the variance of the number of levels in an energy 
interval that contains r levels on average. (Note that 



3 



the unfolded spectra is dimensionless with average level 
spacing equal to unity, so an energy interval of length r 
contains r levels on average.) Define n(e,r) to be the 
number of levels in a small energy interval of length r at 
energy e (i.e., the number of levels between e and e + r). 
Then T, 2 (r) is given by 

2 

S 2 (r) = n 2 (e,r) - n(s, r) — n 2 (e,r) - r 2 , (7) 

where the overbar denot es a ru nning average over the 
measured spectra; note n(e, r) = r only for unfolded 
data. 

In most nuclear physics applications, a different statis- 
tic known as the spectral rigidity, A 3 , is preferred over 
£ 2 (r), primarily because A3 is smoother and fits data 
to a GOE more easily [16| . Qualitatively, A3 measures 
a spectrum's deviation from the best-fit uniform spec- 
trum. The two statistics are related however, and in fact 
A3 can be expressed as an integral over £ 2 (r) [l9j]. How- 
ever, A3 also "washes out" the distinction between GOE 
and GUE at small spacings where TOPE would be most 
evident, so for our purposes S 2 (r) is more useful. 

The crucial point is that TOPE interactions make 
S 2 (r) smaller because of greater level repulsion. Lengthy 
derivations [2(| demonstrate that for the relevant case of 
small A, 

S 2 (r,A) «S 2 (r,0)-4A (8) 

i.e., the variance decreases and spectral uniformity in- 
creases as level repulsion is increased from the T- 
preserving case £ 2 (r, 0). 

III. SENSITIVITY AND EXPERIMENTAL 
REALIZATION 

We now consider how stringent a bound such a test 
could place on a. The sample error a is simply a func- 
tion of the number of levels considered p, so from the cen- 
tral limit theorem we would expect to have a oc p~ x / 2 . 
Detailed calculations verify this 20]; in fact, S 2 (l) has 
a Xp distribution, so the sample error is simply a ~ 
0.6(2/p) 1 / 2 S 2 (l,A) 33]. Roughly, E 2 (1,A) « 1/2 [H. 
Therefore, to obtain a bound on a at the 10~ 3 level, com- 
parable to the experimental bound from Eq. and an 
order of magnitude more stringent than the positronium 
bound from Eq. (QJ, would require ~ 10 5 levels. 

How many chaotic atomic states can we expect to find? 
Obviously there are infinitely many levels of a given J and 
parity including Rydberg series, but we are only inter- 
ested in chaotically mixed compound states. The authors 
of Ref. [2l| performed configuration interaction calcula- 
tions on Th in order to estimate density of states. For Th 
II, they estimate that there are ~ 10 3 compound states 
with J = 3/2 and even parity. Except for a few low- lying 
states, nearly all of these levels are chaotic and suitable 
for our purposes. We assume a comparable number of 
levels would be found for different angular symmetries 



and in different atoms. Spectra could be measured for 
the neutral and singly ionized lanthanides (perhaps even 
doubly- ionized) , giving perhaps 10 — 30 species depend- 
ing on experimental challenges. If ~ 10 sequences with 
different angular symmetry are measured in each species, 
achieving > 10 5 levels seems optimistic but plausible. 

Accurately locating ~ 10 5 levels presents a daunting 
task, but it should be possible with reasonable efficiency 
using frequency comb (FC) spectroscopy. Prior works 
have used the molecular fingerprinting technique to mea- 
sure molecular spectra in a massively parallel fashion 
(see, e.g., [22H24I ] . or [25|, 26] for surveys of experimental 
methods). Essentially, these experiments use the many 
teeth of the FC like thousands of cw-lasers simultane- 
ously probing the sample. The same technique could be 
applied to the dense level structure of the lanthanides. 
However, for our purposes we need spectra with no spu- 
rious or missing levels, which demands FC coverage from 
the visible to the mid-UV, and perhaps even into the 
vacuum-UV or extreme-UV for ionized species. Such 
FCs already exist in the visible and near-IR, and rapid 
progress has been made in extending frequency combs 
to the mid-IR [13], to the VUV [H], and even into the 
XUV [29M3l1 ] . As this technology continues to improve, 
measuring the necessary spectra seems quite achievable. 

Even with FC spectroscopy, making accurate angular 
momentum and parity assignments for so many levels 
would be challenging. It has been suggested [20| that a 
Bayesian analysis could help in this task, since the ap- 
proximate form for the spacing distribution is known. 
This would introduce a few spurious and missing lev- 
els. Depending on the experimental challenges, however, 
it may be useful for labeling 10 5 levels. Studies in nu- 
clear spectra have shown that the effect of spurious and 
missing levels can be accounted for, and as long as their 
number is sufficiently small, they do not invalidate the 
predictions of RME theory [ill [20| ■ 

Finally, we note that this proposal is best suited to a 
search for nuclear-spin independent TOPE interactions, 
experimentally and theoretically. Experimentally, the 
bound on nuclear-spin dependent TOPE from Eq. ([2]) 
is currently an order of magnitude stronger than that on 
e-e TOPE interactions from Eq. (p}. The theoretical rea- 
sons are as follows. When we consider the Hamiltonian to 
be block-diagonal, there is some ambiguity as to what we 
consider to be the total angular momentum. First, let us 
neglect the nuclear spin and the hyperfine interaction. In 
the lanthanides, total orbital and spin angular momenta 
L and S are not "good" quantum numbers thanks to 
strong spin-orbit coupling. In fact, the total electronic 
angular momentum J = L + S is the only good quan- 
tum number, along with parity, so we may assign blocks 
according to a particular J and parity. Nearest-neighbor- 
spacings within a block will follow the Wigner distribu- 
tion, indicating chaos [li| . One may then compute 
fluctuation measures (i.e., £ 2 (r)) in each block, searching 
for deviations from GOE statistics. Any detected TOPE 
interaction would clearly be nuclear-spin independent. 
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Alternately, if we include hyperfme interactions, the 
true total angular momentum of the system is F = J + 1, 
where / is the nuclear spin. Then we may choose blocks 
according to a particular F and parity. However, J is still 
very nearly a good quantum number, which qualitatively 
means that states of a particular F are not chaotically 
mixed. As shown in Ref. [l8j], this means that the near- 
est neighbor distribution for a particular F block will not 
be the Wigner distribution, but rather a superposition of 
several Wigner distributions corresponding to the differ- 
ent possible J values. This distribution will interpolate 
between a Wigner and exponential distribution. Since 
the level statistics would not correspond to the GOE 
even approximately, it would be impossible to search for 
a TOPE interaction represented by a small admixture 
of GUE. While the level positions would certainly be af- 
fected by the existence of some nuclear-spin-dependent 
TOPE interaction, such an interaction would not be de- 
tectable from level statistics using the methods of this 
paper. 

We would like to thank V. Dzuba, V. Flambaum, G. 
Gribakin, M.G. Kozlov, J. Lawler, K. A. Mitchell, and 
M. Pospelov for discussions. This work was supported in 
part by the NSF. 

Appendix 

In this appendix we give an illustrative derivation of 
the behavior of the level-spacing distribution for the T- 



even and T-odd Hamiltonians for small level separations. 
As long as we are interested in very small energy split- 
tings A between neighboring levels, we can neglect all 
other states and consider just these two levels. The 
random-matrix Hamiltonian then reads, 



For the T-even Hamiltonian the parameter c can be set 
to zero, while for the T-odd Hamiltonian all three pa- 
rameters are nonzero. By diagonalizing this matrix we 
find the splitting to be 



A = 2\/a 2 + & 2 + c 2 . 



Let us assume that all nonzero parameters are ran- 
dom variables with zero average and the same variance 
a 2 . Then the probability that A < e <C o is ~ (e/c) 3 
for T-odd Hamiltonian and ~ (e/cr) 2 for T-even Hamil- 
tonian. The probability density p(e) is then ~ (e/cr) 2 
and ~ (e/cr) respectively. Clearly, the level statistics is 
affected by underlying symmetry of the Hamiltonian. 
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